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Electrical microstimulation studies provide some of the most direct evidence for the 
neural representation of muscle synergies. These synergies, i.e., coordinated activations 
of groups of muscles, have been proposed as building blocks for the construction of motor 
behaviors by the nervous system. Intraspinal or intracortical microstimulation (ICMS) 
has been shown to evoke muscle patterns that can be resolved into a small set of 
synergies similar to those seen in natural behavior. However, questions remain about the 
validity of microstimulation as a probe of neural function, particularly given the relatively 
long trains of supratheshold stimuli used in these studies. Here, we examined whether 
muscle synergies evoked during ICMS in two rhesus macaques were similarly encoded 
by nearby motor cortical units during a purely voluntary behavior involving object reach, 
grasp, and carry movements. At each microstimulation site we identified the synergy 
most strongly evoked among those extracted from muscle patterns evoked over all 
microstimulation sites. For each cortical unit recorded at the same microstimulation site, 
we then identified the synergy most strongly encoded among those extracted from 
muscle patterns recorded during the voluntary behavior. We found that the synergy most 
strongly evoked at an ICMS site matched the synergy most strongly encoded by proximal 
units more often than expected by chance. These results suggest a common neural 
substrate for microstimulation-evoked motor responses and for the generation of muscle 
patterns during natural behaviors. 
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INTRODUCTION 

In numerous studies, motor primitives have been defined as "syn- 
ergies" in which a group of muscles is simuhaneously recruited, 
each with a specific balance of activation. These investigations 
have involved frog axial and hindlimb behaviors (Tresch et al., 
1999), cat axial and hindlimb behaviors (Ting and Macpherson, 
2005) and forelimb reaches (Yakovenko et al, 2010), rat fore- 
limb reaches and grasps (Kargo and Nitz, 2003), monkey forelimb 
reaches and grasps (Brochier et al, 2004; Overduin et al, 2012), 
and human axial and hindlimb behaviors (Torres-Oviedo and 
Ting, 2007) and forelimb reaches, grasps and gestures (d'Avella 
et al, 2006; Klein Breteler et al, 2007; Berger et al, 2013). 
Analogous synergies have also been defined at the kinematic level 
(e.g., Santello et al, 1998; Centner and Classen, 2006; Centner 
etal, 2010). 

Much of the direct evidence for synergistic muscle control 
by the central nervous system (CNS) comes from studies of the 
spinal cord. For instance, chemical intraspinal microstimulation 
in the frog evokes topographically-organized, low-dimensional 
electromyographic (EMC) activity patterns (Saltiel et al., 2001). 
In spinalized frogs, too, neurons including those in the interme- 
diate zone of the spinal cord are better correlated with synergistic 
premotor drives than to the activity of individual muscles (Hart 



and Ciszter, 2010). The spinal cord of primates contains premotor 
interneurons facilitating multiple muscles, including those intrin- 
sic to the hand (Takei and Seki, 2010), and has been proposed 
as a substrate for synergies (Tresch et al, 1999; Cheung et al., 
2009). Phasic activation of such units may be responsible for 
multi-muscular EMC bursts (Kargo and Ciszter, 2008). 

Synergies may be encoded at supraspinal as weU as spinal lev- 
els, e.g., in the brainstem (Roh et al., 201 1). At the level of the pri- 
mary motor cortex (MI), in rodents learning a reaching task the 
firing rates of a minority of neurons are correlated with changes in 
the activation of synergies extracted from forelimb muscle EMC 
activity (Kargo and Nitz, 2003). Continuous neural control of 
synergistic muscle groups in primates is also circumstantially sug- 
gested by the ability to reconstruct foreUmb EMC profiles as the 
weighted sum of an ensemble of neurons (Morrow and Miller, 
2003; Schieber and Rivlis, 2007). MI neurons have muscle fields 
(defined by the strongest cell-EMC correlations) that appear to 
fall into relatively few, synergy-like clusters (Holdefer and Miller, 
2002). Such cortical muscle fields may be "hard-wired," chang- 
ing their structure rarely, if ever (Kargo and Nitz, 2003). This 
may be particularly true of corticospinal and corticomotoneu- 
ronal cells, which facilitate small sets of muscles relatively directly 
(Fetz and Cheney, 1980; Bennett and Lemon, 1994). The latter cell 
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population may define a "new MI" that affords muscular coor- 
dination unconstrained by synergies encoded in the spinal cord 
(Rathelot and Strick, 2009). 

As with intraspinal microstimulation, application of relatively 
long trains of suprathreshold electrical current to individual 
sites in motor cortex can evoke complex movements. In cats 
(Ward, 1938), rats (Haiss and Schwarz, 2005; Ramanathan et al., 
2006), prosimians (Stepniewska et al, 2005, 2011), and macaques 
(Graziano et al., 2002a, 2004b, 2005), suprathreshold microstim- 
ulation trains lasting several hundred milliseconds evoke complex 
multijoint forces that frequently drive the animal's body toward 
invariant postures. Shorter- train (<100 ms) intracortical micros- 
timulation (ICMS), in contrast, typically evokes simpler twitch- 
like movements, often restricted to single joints (Graziano et al, 
2002a; Stepniewska et al., 201 1). 

In the case of forelimb and hand areas of non-human primate 
motor cortex, ICMS has been shown to evoke behavioral frag- 
ments including reaching and defensive motions (Graziano et al., 
2002a; Kaas et al, 2013). It has been suggested that motor cortical 
areas may be defined by a continuous map of endpoint postural 
space (Graziano et al, 2002b), e.g., divided according to the 3D 
regions around the monkey to which the forelimb is driven by 
ICMS (Graziano et al, 2002a), or to distinct behaviors like reach- 
ing and defending (Graziano et al., 2005). Whatever the nature 
of the topographical clustering of ethologicaUy-relevant move- 
ments on the motor cortical surface, this organization appears 
to be reflected by distinct, interconnected regions in premotor 
and posterior parietal cortex (Stepniewska et al., 2005, 201 1, 2013; 
Gharbawie et al., 2011). 

Recently, we demonstrated that ICMS of monkey motor cortex 
elicited EMG patterns that could be decomposed into mus- 
cle synergies — ones similar to those seen in natural behavior 
(Overduin et al., 2012; for discussion see also Diedrichsen and 
Classen, 2012; Bizzi and Cheung, 2013; Santello et al, 2013). 
These EMG patterns co-occurred with movements of the hand 
toward a convergent posture. Like the apparent topographical 
organization of spinaUy-encoded synergistic limb movements 
(Tresch et al, 1999; Bizzi et al, 2000, 2002, 2008; Tresch et al, 
2002), we also observed a non-uniform representation of each 
forelimb synergy on the cortical surface of macaques (Overduin 
et al., 2012). Much as microstimulation of multiple points in the 
spinal cord of frogs (Mussa-Ivaldi et al., 1994; Lemay et al., 2001) 
and rats (Tresch and Bizzi, 1999) evokes a linear summation of 
convergent forces and end posture, when ICMS is applied at mul- 
tiple points in the motor cortex of anesthetized cats, the evoked 
EMG activity tends to sum linearly (Ethier et al., 2006). Together 
these results suggest a role of muscle synergies in simplifying 
(even to the point of linearizing) motor control. 

Yet even if ICMS is able to evoke movements through com- 
bination of hard-wired muscle synergies, this does not imply 
that the synergies are organized intracortically. In humans, for 
instance, cortical stroke appears to spare some or most synergies 
(Cheung et al, 2009, 2012; Cruz and Dhaher, 2009), suggest- 
ing encoding at a subcortical locus. In the absence of descending 
signals from the brain, the spinal cord remains fuUy capable of 
generating complex behaviors and muscle activations (Pearson 
and Rossignol, 1991; Zimmermann et al, 2011), as well as 



microstimulation-evoked convergent forces (Giszter et al., 1993; 
Aoyagi et al, 2004). 

Here, we examined motor cortical activity during voluntary 
behavior to see if it might play any role in controlling the 
activation (if not the structure) of muscle synergies. In par- 
ticular, we tested a simple experimental hypothesis motivated 
by our earlier microstimulation work (Overduin et al., 2012), 
namely that cortical units should preferentially encode syner- 
gies similar to those evoked by ICMS at the same electrode. 
(Our null hypothesis, in contrast, was that whichever synergy 
was most strongly encoded by a given unit would bear only 
at-chance similarity to those evoked by ICMS near the unit.) 
If so, this would suggest that intracortical currents — whether 
endogenously generated by motor planning or exogenously intro- 
duced by microstimulation — may determine the degree to which 
downstream synergies are recruited. 

MATERIALS AND METHODS 
SUBJECTS 

Behavioral, muscular, and cortical data were collected from two 
rhesus macaques {Macaca mulatto): Gl (a 5.9-kg, 8-year-old 
female) and G2 (a 6.5-kg, 4-year-old male). Procedures were 
approved by the MIT Committee on Animal Care, and conformed 
to the National Institutes of Health Guide for the Care and Use of 
Laboratory Animals. 

BEHAVIOR 

Subjects used their left hand to press a start button and then 
reach for, grasp, and carry objects between two wells spaced 20 cm 
apart. Reward (0.2-0.3 ml of apple juice or water) was given if the 
object was removed from the first well within 1.0 s, and released 
into the second well within another 1.0 s, where it had to remain 
for at least 0.1 s. Button press, reward dispensation, and data from 
two photosensors (E3T-SR12; Omron, Kyoto, Japan) mounted 
within each well were recorded together, allowing trials to be 
divided into reach and carry phases (Overduin et al, 2008). The 
25 Delrin plastic objects (density 1.4g/cm^) included 5 spheres 
of variable diameter (ranging from 1.6 to 3.6 cm), 5 cubes of 
variable width (1.5-3.6 cm), and 15 cylinders of which 5 each 
spanned one of three dimensions (height, 0.6-5.7 cm; uniform 
diameter, 1.3-3.8 cm; inner diameter, 0.6-3.2 cm). The number 
of objects presented in a given day, or spanned by any one single 
unit recorded in a given session, could be fewer than 25. The same 
object was presented enough times consecutively for the animal 
to be able to perform 10 successful trials, before another object 
was pseudorandomly selected. During recordings, subjects were 
head-restrained via an implanted cranial post. 

SESSIONS 

Cortical recordings from monkey Gl comprised 7798 successful 
left-target-directed trials performed over 20 recording sessions 
spanning 45 days; those from G2 comprised 775 left-target- 
directed trials performed over 6 sessions spanning 6 days. The 
analysis presented here, however, was focused on the subset of 
these data for which ICMS experiments were performed on the 
same days as the cortical recordings (4485 trials over 13 ses- 
sions from Gl and 544 trials over 4 sessions from G2). Muscle 
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recording was done in each of these sessions from Gl and G2, 
and also included trials (as well as other, interspersed sessions) 
when cortical recording was not done. The EMG data from which 
synergies were extracted (Overduin et al., 2008, 2012; Figure 3A) 
thus include trials recorded both with and without simultaneous 
neural data. In particular, of Gl and G2's 1000 left-target-directed 
trials used to construct trial-averaged EMG activity for synergy 
extraction, 482 and 462 trials overlapped with the 4485 and 544 
trial subset we focus on here, respectively. The EMG data of the 
remaining 89% (4003/4485, Gl) and 15% (82/544, G2) of trials 
are new to this report. None of the cortical data have previously 
been presented. 

SURGERY 

Cortical surgeries followed (Gl) or occurred along with the first 
of (G2) the muscle electrode implantation surgeries (described 
in Overduin et al., 2008). These surgeries were performed under 
sterile conditions and general anesthesia (0.05 mg/kg atropine 
and 10 mg/kg ketamine injected intramuscularly, followed in Gl 
by 5 mg/kg sodium pentobarbital intravenously or in G2 by 
inhalation of 1-2% isoflurane with 2 L O2). Craniotomies were 
centered over right-hemisphere motor cortex. Custom stainless 
steel wells (Gl: 28 mm wide, G2: 20 mm) were secured with 
bone screws and bone cement. The animals were given analgesics 
and systemic antibiotics following the surgeries. The dura was 
kept intact during surgery, and was subsequently treated with 
topical antibiotics and anti-inflammatories. Fresh connective tis- 
sue growth above the dura was further controlled by periodic 
(~ Ix weekly) mechanical scraping, done under light anesthesia 
through the weeks of recording. 

MUSCLES 

EMG recordings were made via 15 (Gl) or 19 (G2) elec- 
trodes chronically implanted in muscles of the left forelimb. 
Proximal muscles acting on the shoulder and elbow included 
Del (deltoideus). Pec (pectoralis major), TriU and TriR (tri- 
ceps brachii, ulnar and radial short heads), Bic (biceps brachii 
longus), and BR (brachioradialis). Wrist and extrinsic hand 
extensors included AbPL (abductor poUicis longus), ECRB 
(extensor carpi radialis brevis), EDC (extensor digitorum com- 
munis), ED23 (extensor digiti secundi and tertii proprius), ED45 
(extensor digiti quarti and quinti proprius), and ECU (exten- 
sor carpi ulnaris). Wrist and extrinsic hand flexors included 
FCR (flexor carpi radialis), EDS (flexor digitorum superfi- 
cialis), FDPU and FDPR (flexor digitorum profundus, ulnar and 
radial), and FCU (flexor carpi ulnaris). Intrinsic hand mus- 
cles included AbPB (abductor poUicis brevis), AdP (adductor 
pollicis), OpP (opponens poUicis), F5B (flexor digiti quinti bre- 
vis manus), and Op5 (opponens digiti quinti manus). EMG 
data were recorded on a trial-by-trial basis (between button- 
press and reward events). These data were bandpass-filtered 
(between 10 and 1000 Hz), notch filtered (60 Hz), and differ- 
entially amplified (5000x) by a programmable signal condi- 
tioner (CyberAmp 380; Molecular Devices, Sunnyvale, CA) under 
software control (CyberControl; Molecular Devices). Data were 
digitized (2 kHz) via a data acquisition board (NI PCI-6035E; 
National Instruments, Austin, TX) under custom software control 
(LabVIEW; National Instruments). EMG channel subselection 



following cross-talk analysis is described in Overduin et al. 
(2008). 

CORTEX 

Single units were recorded from dorsal premotor cortex (PMd), 
ventral premotor cortex (PMv), and MI (Figure 1). Areas were 
identified using magnetic resonance imaging data and sensori- 
motor mapping including ICMS. MRIs were collected 2 years 
after (Gl) or 4 months before (G2) the recordings presented 
here. At the beginning of recording sessions, unit somatosensory 
(proprioceptive and cutaneous) response fields were estimated by 
passively moving the monkey's limbs and by stimulating its skin. 
Both in early exploratory sessions and at the end of each recording 
session, ICMS was applied via tungsten microelectrodes (FHC, 
Bowdoin, Maine; 0.3-3-MQ impedance, 250-|J-m shaft diame- 
ter tapered to a 3-|im-wide tip). Up to 10 such electrodes were 
acutely introduced into the brain in each session using a man- 
ual microdrive (30-|xm depth resolution). The microdrive was 
mounted on a grid that was secured to the recording well and 
that constrained the interelectrode spacing to 1 mm. The ICMS- 
evoked movements and stimulation thresholds were used to iden- 
tify the portion of cortex sampled by the electrodes. Stimulation 
parameters used for this mapping included 2 x 0.2-ms pulse 
duration (cathodal-leading), 10-150-|iA current, 330-Hz pulse 
frequency, and 0.05-s train duration. Modified ICMS parameters 
were used to evoke longer-lasting movements at the end of some 
sessions (Figure 2A): 8-100-|xA current, 200-Hz pulse frequency, 
and 0.15-0.5-s train duration (as in Overduin et al., 2012). We 
consider only the first seven ICMS trains delivered at each site, 
this being the minimal number applied across the 33 (Gl) or 13 
(G2) sites (including 32 in MI, 9 in PMd, and 5 in PMv; Overduin 
etal, 2012). 

UNITS 

Extracellular, intracortical voltages were recorded through the 
same electrodes used for ICMS. During recordings, no attempt 
was made to record from units within particular cortical layers 
or in sulcal sites, and the laminar location of recorded units was 
generally unknown given limited depth resolution and dimpling 
of the cortex upon penetration. Instead, we selected for record- 
ing the first stably-firing and well-discriminated unit(s) (if any) 
discovered along a given electrode track, as revealed by audio and 
oscilloscope monitoring of the recorded voltage. We rejected all 
cortical sites wherein somatosensory stimulation or ICMS had 
indicated response fields extending beyond the left forelimb (e.g., 
to the face or leg). We did not require the remaining sites to 
demonstrate sensorimotor responses, only that they be within 
topographical regions convexly bounded by sites that did. No 
other criteria were applied to units; i.e., we treated all as poten- 
tially task-related. We only consider units recorded at ICMS elec- 
trode sites (33 sites in Gl, 13 in G2; Overduin et al., 2012); these 
comprise 94 units (83 in Gl, 11 in G2). Signals were preamplified 
at unity gain by a headstage located ~5 cm from the electrodes, 
and then passed to an amplifier for amplification ( 10,000 x ) and 
bandpass-filtering (600-6000 Hz, 2nd-order filter with roll-off on 
both ends) before digitization (Neuralynx, Inc., Tucson, Arizona). 
Spikes were identified online when electrode voltages exceeded 
a manually-set threshold, and 1.1 -ms waveforms (including the 
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threshold-crossing moment at 0.26 ms into the waveform), sam- 
pled at 30 kHz, were stored to disk. Spike times and other 
event times were recorded together for later synchronization of 
behavioral, EMG, and neural data. Offline, single units were 
identified based on spike waveform features and interspike inter- 
vals (ISIs), using manual clustering with MClust (MClust-3.4, 
A. D. Redish et al.) and custom routines written in MATLAB 
(Math Works, Natick, Massachusetts), for this and the following 
analysis. Particular care was taken to ensure that waveform fea- 
tures and firing rates remained relatively constant over the record- 
ing span of each accepted unit. The 94 units had mean firing rates 
between 1.0 and 58.0 Hz, and signal-to-noise ratios between 2.4 
and 24.3 (deCharms et al., 2009), with 1% of ISIs <1 ms. 

PREPROCESSING 

After preprocessing as described in Overduin et al. (2008), 
grasping- related EMG data were integrated over 9-ms (Gl) or 
11-ms (G2) bins, and normalized to each channel's maximum 
integrated EMG over object conditions. The EMG data of 40 
trials within each of the 25 object conditions were aligned on 
the time of object removal from the first well, cropped to a 
100-sample window around this event [Gl:(— 0.35 : -1-0.55) s, 
G2:(— 0.5 : -1-0.6) s], and then averaged over trials within each 
object condition. (The different window and integration times 
for the two animals were chosen based on their different move- 
ment latencies; Overduin et al, 2008.) ICMS-evoked EMG data 
were integrated over a (-1-0.025 : -1-0.150) s window relative to 
ICMS onset for each of 7 trains applied at each stimulation 
site, and normalized by the same factors as the grasping-related 
EMG data (Overduin et al, 2012). With regards to the neural 
data, mean firing rates were computed on a within-object basis, 
using all trials fully spanned by a unit. Trials were time-aligned 
on the time of object removal from the origin well (as for the 
grasping-related EMG data). Subsequent analysis was restricted to 
a (-0.4 : -1-0.5) s (Gl) or (-0.55 : -1-0.55) s (G2) window around 
this event, i.e., to the same windows as for the grasping-related 
EMG data minus a fixed 50-ms delay (Morrow and Miller, 2003; 
Schieber and Rivlis, 2007; Stark et al, 2007). Within this window, 
each unit's spikes were summed within 9- or 11-ms bins (again, 
for consistency with the grasping-related EMG data binning). The 
units' mean firing rate profiles within each object condition were 
then smoothed by convolution with a 50-ms Gaussian kernel. 
(This smoothing, evident in Figure 3B, did not have a qualita- 
tive effect on the summary results of Figure 4.) All preprocessing 
was done in MATLAB. 

SYNERGIES 

Non-negative matrix factorization (NNMF) was used to identify 
a set of synchronous synergies underlying each monkey's mus- 
cle patterns (Lee and Seung, 1999; Tresch et al., 1999). There were 
two sets of muscle patterns to consider: ICMS-evoked EMG data / 
(including the data evoked at the sample site in Figure 2A, right), 
and grasping-related EMG data G (shown averaged over object 
conditions in Figure 3A, fooftom). Matrix I(e, s, Z) pools together 
the activity in E EMG channels evoked by each of S stimulation 
trains delivered at each of L ICMS locations (Overduin et al., 
20 12). 7 thus has dimensionality 15 x 7 x 33(Gl)orl9 x 7 x 13 



(G2). Matrix G(e, f, o) specifies trial-averaged EMG activity 
over the same E channels, but across each of T time points 
in each of the O object conditions. G thus has dimensionality 
15 X 100 X 50 (Gl) or 19 x 100 x 50 (G2). Data I and G were 
independently reconstructed as combinations of n = 1 , • • • , 
or M = 1 , • • • , N'^ synergies, respectively, where both and N'^ 
were evaluated up to the number of EMG channels E, i.e., 15 (Gl) 
or 19 (G2). Each ICMS-evoked or grasp-related synergy vj,(e) 
or v^(e) is a vector of length E capturing a unique balance of 
activation across the EMG channels. In reconstructions, each syn- 
ergy was multiplied by a scalar, non-negative weighting coefficient 
w|j(s, I) or wj^(f, o). These coefficients could vary both within 
conditions (i.e., over stimulation trains s or time samples f ) and 
across conditions (i.e., over ICMS locations / or objects o). The 
reconstructions can be expressed as: 

I(e, 5,0= ^wi„(5, /).v!,(e) (1) 

n= 1 

G{e, t, o)=J2^n(t, o)-v'^(e) (2) 

n= 1 

For a given or N'^, the algorithms iteratively updated the struc- 
tures vj, or and coefficients w), and wj^ until the total recon- 
struction error, R^, increased by less than 0.001 over 10 iterations. 
The algorithms were repeated five times for each extraction; the 
set of synergies with the highest EMG variation explained was 
selected for further analysis. Dimensionalities (N^ and N'^) were 
chosen by applying a threshold of = 95% (Overduin et al., 
2012). 

MATCHING 

Synergies vj, were compared and matched to synergies using 
a greedy search procedure (Tresch et al, 1999; Overduin et al., 
2012). For all JV' x N'^ possible pairs of ICMS-evoked and 
grasping-related synergies, we first computed dot products (e.g., 
6 X 8 = 48 dot products, for G2. The pair of ICMS-evoked 
vs. grasping-related synergies with the highest dot product was 
defined as the best-matching pair. The pair with the highest dot 
product among the remaining (N^ — 1) x (N'^ — 1) pairs (e.g., 
5 X 7 = 35 pairs, for G2) defined the second-best match. This 
process was repeated until all synergies in one set had been paired 
(e.g., over min(6,8) = 6 times for G2). We then used Monte Carlo 
simulation to assess the significance of each match. We repeated 
the greedy search algorithm 10,000 times for each monkey, after 
first randomly shuffling EMG channel identity each time. (For 
G2, for instance, this involved finding 10,000 x 48 = 480,000 
dot products.) We then compared the highest (best-matching) 
dot product between actual ICMS-evoked and grasping-related 
synergies with the distribution of highest dot products from the 
10,000 comparisons of shuffled synergies. If the former value 
exceeded the 95th percentile of the distribution of latter values, 
we took the match as significant at p < 0.05. We then repeated 
this comparison for the second-best actual synergy pair vs. the 
distribution of second-best shuffled pairs, etc. 
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FIGURE 1 I A portion of each monkey's right hemisphere as viewed 
from above. Solid lines represent sulci (CS, central sulcus; SPcS, 
superior precentral sulcus; ArS, arcuate sulcus); daslied lines show 
estimated inter-areal borders separating Ml, PMd, and PMv. The circles 
show stimulation and recording sites, with dots indicating the number 
of units (if any) recorded at the site. For G2, ICMS sites are colored 
according to the degree to which the evoked muscle EMG activity was 
similar to (i.e., reconstructed by) synergy vj,. (Specifically, the color 
scale indicates the mean amplitude coefficient used in 
reconstructing the ICMS-evoked EMG vectors at a site with synergy yj,, 
averaged over ICMS trains.) Of G2's ICMS sites, this synergy was 
most strongly evoked at site XI, in Ml. 



ANALYSIS 

In reconstructing the ICMS-evoked EMG vectors at a given site, 
the synergy vj, most "evoked" at the site was the one with 
the largest weighting coefficient wj,, averaged over ICMS trains 
(Figure 2A, showing results for site XI in Figure 1, right). To 
determine which synergy was instead most "encoded" by a unit, 
the unit's firing rate profile was correlated against the task-related 
synergy scaling profiles w'^, and the largest positive correla- 
tion was identified. Results were insensitive to the use of linear 
(Pearson) or rank (Kendall or Spearman) correlation. To deter- 
mine which task-related synergy was most evoked at a recording 
site, the site's ICMS-evoked EMG patterns were decomposed into 
combinations of synergies, which in turn were matched to task- 
related synergies (as described above). Over units, we counted the 
frequency at which the most-encoded synergy was the same as 
the most-evoked synergy at the electrode, defined as above. The 
chance frequency of such matches was 1/6 = 17%, as the analy- 
sis considered the 6 muscle synergies both evoked by ICMS (in 
the case of G2, Figure 2A, left) and observed in the task data 
(Figure 2B). Results of this x.^ test were evaluated for significance 
at ap < 0.05 threshold. 

RESULTS 

In a recent study (Overduin et al., 2012), we examined ICMS- 
evoked movements and muscle activity in two rhesus macaques 
("Gl" and "G2"). EMG data were recorded from 15-19 electrodes 
chronically implanted in muscles of the shoulder, arm, and hand. 
As we reported (Overduin et al., 2012), ICMS within MI and 
dorsal and ventral premotor (PMd and PMd) cortex (Figure 1) 
appeared to drive the forelimb toward an invariant, site-specific 
posture, and at the same time to replace voluntary muscle activity 
with an invariant, site-specific tonic EMG pattern (as in Griffin 
etal, 2011). 

After aggregating ICMS-evoked EMG data over multiple stim- 
ulation sites, these pooled data could be decomposed into 
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FIGURE 2 I Microstimulation-evoked muscle activity and 
reconstruction by synergies. (A) The EMG activity evoked during the first 
0.15s of ICMS at site XI from Figure 1, shown both in gold-colored profile 
{far right) and as the integrated mean ± SE of this activity over stimulation 
trains {second from right). (The vertical 20-|j,V scale bars to the right of the 
EMG profiles give the relative level of activation in each channel.) The mean 
activity could be reconstructed by six ICMS-evoked synergies vj, {left), with 
the largest weighting given to v],. (B) Of eight muscle synergies v'^ derived 
from an object grasping behavior, synergies v^-vg could be matched in 
structure to ICMS-evoked synergies v^-vg. 



combinations of a reduced set synchronous synergies using 
NNMF (Lee and Seung, 1999; Tresch et al, 1999). Each of the 
ICMS-derived synergies v\ (n = l,...,N^) captures a pattern of 
synchronous firing over muscles. We found that N' = 6 (G2) or 
7 (Gl) synergies were sufficient to reconstruct > 95% of the vari- 
ability in the ICMS-evoked EMG data (Overduin et al., 2012). 
Such data reconstruction is exemplified in Figure 2A for G2's 
ICMS site XI, in MI. The figure depicts the ICMS-evoked activ- 
ity at this site in the form of: the time-varying, trial-averaged 
EMG signal (far right); an ICMS-evoked EMG vector integrating 
this per-ICMS activity over time (middle); and as a vector sum 
of ICMS-derived synergies, each scaled by the weighting coeffi- 
cient w\ (left). For this site, the ICMS-evoked EMG activity was 
dominated by synergy v^, as shown. 

We also studied a manual behavior performed by these sub- 
jects (Overduin et al., 2008, 2012), in which they had to reach 
for an object presented in a well, grasp it, and then carry it to 
the opposing well. To elicit a variety of hand postures and forces, 
the objects included 25 spheres, cubes, and cylinders of differ- 
ent dimension. Pooled over multiple days, the EMG data could 
be decomposed into combinations of grasping-related synergies 

(n = 1,. . . ,N^). We found that = 8 (G2) or 10 (Gl) syn- 
ergies were sufficient to reconstruct > 95% of the variability in 
the grasping-related EMG data (Overduin et al., 2012). For each 
monkey, 6 of these synergies could be matched uniquely to one 
of the ICMS-derived synergies (Figures 2A vs. 2B) using a greedy 
search procedure. Figures 2B, 3A illustrate the synergies found 
for one animal and their reconstruction of average EMG activity. 
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FIGURE 3 I Task-related muscle activity and reconstruction by 
synergies, vs. cortical activity. (A) EMG activity in the grasping task (gray 
envelopes at bottom] could be decomposed as a summation of 
grasping-related synergies vj; {red lines at bottom), each weighted by the 
coefficients [top). EMG activity and synergy coefficients are shown as 
mean {bottom) and mean ± SE (top), over all 25 objects and 733 trials 
spanned by this monkey's EMG dataset. (B) The task-related firing rate of a 
unit recorded at monkey G2' site XI, shown both as a spike raster {bottom) 
and as the binned activity mean ± SE (over the 22 objects and 113 trials 
spanned by this unit, fop). Of all the task-related synergies v^-Vg used in 
reconstructing EMG activity unit's activity (A, top), the one corresponding 
to vlj (i.e., vjj") had an activation profile most similar to the unit's activity 
(B, fop). 



Synergies could provide the animal with a mechanism to 
continuously and efficiently control its muscles, by specifying 
a task-specific amplitude coefficient time course ivj^ for each 
of the synergies (Figure 3A top vs. bottom). The involvement 
of motor cortex in this control is circumstantially suggested 
by the correspondence in synergy structure whether these are 
derived from ICMS-evoked or voluntarily-generated EMG data 
(Figures 2A vs. 2B). Here, we sought more direct evidence for the 
role of motor cortex. In particular, we looked for evidence that the 
muscle synergies evoked by ICMS (and matched to those gener- 
ated voluntarily) were also encoded by single motor cortical units 
near the ICMS electrodes. 

A positive example of such a correspondence is shown in 
Figure 3A (fop) and Figure 3B (fop), for monkey G2. The MI 
unit shown in Figure 3B is at site XI, the location of which is 
highlighted in Figure 1 (right), and the ICMS-evoked EMG activ- 
ity of which is shown in Figure 2A. The single unit recorded 
at this site exhibited a strong burst of activity just prior to the 
time when the object was retrieved from its origin well, and was 
largely quiet during the following carry and release movements 
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FIGURE 4 I Motor cortical units are biased toward encoding synergies 
evoked by nearby microstimulation. For both G1 and G2, the muscle 
synergy most strongly evoked by ICMS was also the synergy most likely to 
be best-correlated with the firing of cortical units at the same electrode. 
The unit depicted in Figure 3B, for instance, contributes to G2's 
solid-colored bar 



(Figure 3B). (These modulations were captured in the unit's fir- 
ing rate, derived from spike counts binned at 50 ms and smoothed 
with a 50-ms Gaussian kernel.) Of all the grasping-related syner- 
gies v^-Vg , the one whose activation profile was most positively 
correlated with that of the unit was (Figure 3A, fop), which 
appeared to capture the coactivation of forearm flexor muscles 
(Figure 2B). And of all the ICMS-evoked synergies v\ — Vg, it was 
the matching synergy v]^ that was most dominant in reconstruc- 
tion of the evoked EMG activity. 

How consistently did the muscle synergy most encoded by 
a unit match the synergy most evoked by ICMS at the same 
site (as in the foregoing example)? For this analysis we sim- 
ply counted the number of MI, PMd and PMv units for which 
the grasping-related synergy most strongly encoded by the unit 
matched the primary muscle synergy evoked by ICMS at the elec- 
trode. (This analysis is restricted to those ICMS sites at which 
units were also recorded; these are highlighted in Figure 1.) 
As shown in Figure 4, the actual frequencies were significantly 
higher than the 1/6 = 17% chance frequency [26 of 94 units, or 
28%; X(i) = 6.82, p < 0.01], as was true for Gl alone [22/83 = 
27% of units; Xq) = 4.82, p < 0.05] and supported by a trend 
among G2's smaller population of units [4/11 = 36% of units; 

= 2.56, p = 0.07]. (Note that the "1st most encoded syn- 
ergy" in this plot is not necessarily vf, but instead whichever 
synergy vj^ was most strongly correlated with the unit's fir- 
ing profile, e.g., for the unit shown in Figure 3B. Similarly, 
the "most evoked synergy" is whichever synergy was dom- 
inant in reconstruction of the EMG activity evoked at the 
unit's electrode. In general, v\, v\, and v\ appeared to be most 
commonly evoked by ICMS, being the synergies "evoked" at 
18/83, 28/83, and 20/83 of monkey Gl's units, respectively 
and by 3/ 1 1 , 3/ 1 1, and 4/ 1 1 of G2's units. ) 

While the qualitative pattern for G2 followed very closely 
that of Gl (Figure 4), it likely failed to reach significance due 
to insufficient sampling of units (11 units, vs. 83 for Gl). 
Also, these trends were weakened by inclusion of PMv units. 
Considering only MI and PMd units, the fraction of cases in which 
the grasping-related synergy most strongly encoded by a unit 
matched the synergy most strongly evoked by proximal ICMS was 
significantly higher than 1/6 in G2 [4/10 = 40% of units; X(i) = 



3.27, p < 0.05] as well as Gl [20/73 = 27% of units; = 5.04 
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p < 0.05], and in both animals combined [24/83 = 29%; = 
7.47,p < 0.01]. 

While the synergies v| extracted from ICMS-evoked EMG 
activity did cluster non-uniformly on the cortical surface 
(Overduin et al., 2012), we observed no tendency of the subset 
of sites at which similar synergies were both encoded and evoked 
to be topographically grouped on the cortical surface. 

DISCUSSION 

Our results suggest that direct stimulation of patches of motor 
cortex generates synergistic combinations of muscle activity that 
are weighted toward synergies encoded at the stimulating elec- 
trode. A common motor cortical substrate appears to be activated 
both by endogenous currents during voluntary behaviors and by 
exogenous currents introduced by electrical microstimulation. In 
either case, these currents appear to modulate the amplitude, 
not structure, of a shared set of downstream muscle synergies. 
Our earlier work also suggests that each synergy (and the pos- 
ture reached through its tonic activation) may be represented by 
a non-uniform map over the motor cortical surface (Overduin 
etal, 2012). 

It may be contested that electrical currents injected into the 
CNS via microstimulation elicit non-natural patterns of neural 
activity (even if the downstream muscle activity is resolvable into 
well-organized movements and natural muscle synergies). The 
biological basis of ICMS, and even the extent of cortex acti- 
vated, remain poorly understood (Butovas and Schwarz, 2003; 
Tolias et al., 2005), despite recent studies using single-cell record- 
ing, behavioral methods, functional magnetic resonance imaging 
(Tehovnik et al., 2006), and two-photon calcium optical imaging 
(Histed et al, 2009). Moreover, these investigations of ICMS have 
largely been limited to more conventional, short-train and/or 
subthreshold ICMS rather than the form used here. In applying 
ICMS in motor cortex, researchers have traditionally used short 
train durations (typically 25-70 ms) and sub- or per-threshold 
currents (typically 10-60 |xA; e.g., Asanuma and Rosen, 1972; 
Sato and Tanji, 1989; Donoghue et al., 1992) in order to map the 
overt response of cell populations including corticomotoneuronal 
cells. Such studies neither sought nor reported convergent move- 
ment responses of the sort found here and by others (Graziano 
et al., 2002a,b), and have naturally aroused some controversy 
(Strick, 2002). 

Together with Overduin et al. (2012), our findings indicate that 
movements, muscle patterns, and cortical activations evoked by 
relatively long-train ICMS can be related to those observed in 
natural behavior. Indeed, it can be argued that a physiologically 
realistic model of motor activation requires ICMS with relatively 
long trains, as well as suprathreshold currents and intermediate 
pulse frequencies (optimally 80-140 |xA and 80-140 Hz, in the 
case of primate MI forelimb-area ICMS; Van Acker et al, 2013). 
Graziano and coworkers (Graziano et al., 2002a) have emphasized 
that ~500-ms stimulation trains approximate the time scale of 
natural movements like primate reach and grasp (Georgopoulos 
et al, 1986; Reina et al., 2001). Even in spinalized amphibians, in 
whom convergent force patterns were first observed (Mussa-Ivaldi 
et al, 1990; Giszter et al, 1993; Loeb et al, 1993), these move- 
ments were evoked by similarly long-train (typically 300-ms) 



intraspinal microstimulation. Researchers using ICMS to study 
oculomotor and somatosensory systems have also used relatively 
long trains. For instance, 400-ms trains were applied to the arcu- 
ate sulcus to replicate the time scale of typical regular head 
movements (Freedman et al, 1996), and 500-ms trains were 
applied to primary somatosensory cortex to mimic a tactile stim- 
ulus (Romo et al, 1998). While trains shorter than ~500 ms elicit 
truncated movements when delivered to primate MI (Graziano 
et al, 2002a; Van Acker et al, 2013), in Overduin et al. (2012) 
we show that even 150 ms of ICMS yields convergent movements 
with a predictable equilibrium point and synergistic muscle 
activation. 

For the majority of units we sampled, firing rate profiles 
were more correlated with muscle synergies other than the most- 
evoked synergy (Figure 4, non-shaded bars). There are many 
ways to account for these other units. For example, units may 
encode other continuously-controlled quantities like the position, 
velocity, or acceleration of extrinsic effectors or intrinsic joints, or 
the forces underlying these kinematics, or the muscle contractions 
determining these dynamics (Carmena et al, 2003; Morrow and 
Miller, 2003; Schieber and Rivlis, 2007; Stark et al, 2007; Velliste 
et al., 2008; Ganguly and Carmena, 2009; Hochberg et al., 2012; 
CoUinger et al., 2013). Motor cortical encoding of purely extrin- 
sic variables may be unlikely based both on theoretical grounds 
(Todorov, 2000; Paninski et al., 2004) and on empirical evidence 
that neurons' intrinsic muscle fields are more stable than their 
tuning to extrinsic dimensions like hand direction (Morrow et al., 
2007), with such tuning observed to fluctuate even during sta- 
ble within-session behavior (Carmena et al., 2005; Rokni et al., 
2007; c£ Chestek et al., 2007). However, the purpose of this report 
was not to determine whether motor cortical units encode muscle 
activity, synergy coefficients, or any other motor variable. Instead, 
we simply sought to test whether the units were more likely than 
chance to encode those synergies evoked by ICMS at the same 
electrode. 

Another way to account for units which appeared to encode 
different synergies than those evoked by nearby ICMS is to 
recognize that the number of units affected by ICMS cur- 
rent far exceeds the handful sampled at an electrode. With the 
stimulation currents needed to elicit complex movements (< 
100 |x A in this report), cortex is no doubt activated far outside 
the 100-|xm radius assumed for ~ IOjjlA stimulation (Ranck, 
1981). Such suprathreshold trains are likely to recruit other 
cells beyond those immediately next to the stimulating elec- 
trode through synaptic connections (Histed et al, 2009), and 
indeed may be required to transsynaptically activate non-direct 
connections between motor cortex and the spinal cord (Strick, 
2002). 

Another possibility to be explored in future work is that pop- 
ulations of units may specify not only continuous variables like 
synergy amplitude, but also sequential recruitment of muscles 
(e.g., via synergies) in the form of discrete "motor programs" 
(Keele, 1968; Polit and Bizzi, 1978, 1979; Georgopoulos et al, 
1983). Elsewhere these programs have been referred to as "time- 
varying synergies," in distinction from the "synchronous syner- 
gies" discussed in the present report (Kargo and Nitz, 2003), 
and have been extracted from muscle data using a modified 



Frontiers in Computational Neuroscience 



www.frontiersin.org 



March 2014 | Volume 8 | Article 20 | 7 



Overduin et al 



Evoked and encoded muscle synergies 



form of NNMF (d'Avella et al, 2003, 2006; Overduin et al, 
2008). A dynamical systems approach may also be appropriate for 
capturing time-varying patterns within population-level neural 
and muscular data (Churchland et al., 2012). 

Besides studies of voluntary neural and muscular dynamics, 
microstimulation studies are also broadly consistent with the idea 
of discrete movement encoding. Intraspinal microstimulation of 
sufficient duration generates forces that tend to drive a limb 
to particular postures or through sequences of postures, in 
frogs (Mussa-Ivaldi et al., 1990; Giszter et al., 1993; Kargo and 
Giszter, 2000), cats (Lemay and Grill, 2004) and rats (Tresch 
and Bizzi, 1999). Sufficiently-long ICMS trains applied to mam- 
malian motor cortex (Ward, 1938; Graziano et al, 2002a, 2004a, 
2005; Stepniewska et al, 2005, 2011; Ramanathan et al, 2006) 
can evoke complex multijoint forelimb behaviors with multiple 
phases (such as reach, grasp, then retraction) and invariant end- 
points. The stimulation-evoked convergent forces (Giszter et al., 
1993; Kargo and Giszter, 2000), invariant endpoints (Graziano 
et al, 2004a; Graziano and Aflalo, 2007), and bell-shaped speed 
profiles (Graziano et al., 2005) all tend to overlap with motions 
and postures found in subjects' natural behavior. Together with 
the present results, these studies indicate that microstimulation- 
evoked movements are not artifactual, and indeed can provide 
insights into natural movement planning. 
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